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A Pegasus computer has been programmed to give the absorption coefficients and virtual heights 
at vertical incidence for the main magneto-ionic components reflected from a stratified ionosphere, 
using the phase integral method. Typical absorption and virtual height curves are given for daytime 
ionospheric models. For high geomagnetic latitudes the Z-trace is shown to be one of the two least 
attenuated traces for daytime flayer reflections at frequencies both above and below the gyrofre- 
quency, and it is shown that under certain circumstances the Z-trace may appear at a lower virtual 
height than the 0-trace. The limitations of the usual ray method which uses the group refractive index 
for determining virtual heights are discussed, and it is shown that a simple correction term may be 
added to the ray theory calculation of virtual heights to make it agree with the phase integral 
determination. 



1. Introduction 

The reflecting properties of the ionosphere have 
been calculated numerically by two main techniques. 
The most accurate is probably the 'full-wave' method 
in which the first order equations derived from Max- 
well's equations and the ionospheric constitutive 
relations are solved by a step-by-step integration. 
The method, described by Budden [1955], Barron and 
Budden [1959] and others, obtains the solutions in 
terms of the electric and magnetic wavefields in per- 
pendicular directions. A recent variant of this tech- 
nique given by Pitteway [1964] considerably simplifies 
the integration procedure by obtaining solutions in 
terms of two independent solutions of the wave equa- 
tion which are not coupled and which in many impor- 
tant cases are similar to the characteristic magneto- 
ionic modes. The full-wave method is suited 
particularly to low-frequency analysis, whereas at 
higher frequencies the large number of steps required 
makes the integration prohibitively long. The results 
obtained by the full wave technique are not readily 
applicable to the propagation of pulses, especially to 
the calculation of virtual heights, since the perpen- 
dicular components in the Budden method, or the non- 
coupled modes in the Pitteway method, must first be 
'unravelled' and reconstituted into characteristic 
magneto-ionic modes in order to permit the differen- 
tiation of phase with respect to frequency which yields 
the virtual heights of the O- and Z-traces. The diffi- 
culty is enhanced when two modes having the same 
polarization, e.g., the and Z modes, exist 
simultaneously. 

The second important technique, the phase integral 
method, has been described by Budden [1961], 
Cooper [1961], and others. It involves integrating 
the complex refractive index n, from ground level 



usually, up to the "height" (usually a complex quan- 
tity) where n becomes zero; i.e., up to one of the 
reflection branch points of the complex refractive 
index. There will generally be two or three such 
reflection points, depending on whether the frequency 
is below or above the gyrofrequency. In addition 
there will be a coupling branch point at a complex 
height where the refractive indices for the character- 
istic magneto-ionic modes become equal. A contour 
of integration which encompasses the coupling point 
(in the sense that the coupling point lies between the 
contour and the real axis) will represent a coupling 
transition from one characteristic mode of polariza- 
tion to another. Not all paths of integration connect- 
ing the reflection points to ground are physically 
significant [Budden, 1961]. In general there will 
be four distinct contours of integration associated 
with modes of reflection which will give independent 
traces on an ionogram — provided the absorption losses 
are not too high. The virtual height of reflection is 
deduced by evaluating the real part of Jndz over these 
integration paths for two adjacent frequencies, and 
the imaginary part of fndz gives the overall absorption. 
The solutions thus give the type of information which 
is directly comparable with the results of pulse re- 
flection experiments and may be used to compare 
the properties of theoretical models of the ionosphere 
with those actually observed. 

2. Theoretical Considerations 

2.1. Appleton-Hartree Formula and the Branch 
Points 

The complex refractive index n is given by the 
Appleton-Hartree formula, which is used in the form 



n 2 =l~- 



X(l-X-iZ) 



1 On leave of absence from the Technion, Israel Institute of Technology, Haifa, Israel. 



{\-iZ)(\-X-iZ)-Y 2 



sin 2 9 



Vf 



(2.1) 
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/here 



C, at a complex height where 2 



Y 4 sin 4 
F a + Y 2 cos 2 0(l-X-iZ) 2 . (2.2) 



The symbols are the usual magneto-ionic parameters 

Ne 2 



X=l + i(Z c -Z). 



(2.6) 



X = 



€ m(27Tf) 2 



in rationalized units, e and 



m being the electronic charge and mass respectively, 
/ the wave frequency, €o the electric permittivity of 

free space, and 
N the number of electrons per unit volume. 

Y =—p= where fn is the gyrofrequency and B 

J rn 

is the induction of the earth's magnetic field. 

Z =- — - where v is the electronic collision frequency. 
Zttj 

6 is the inclination of the earth's magnetic field to 
the vertical. 

In formula (2.1) the square root of the complex quan- 
tity F is chosen using the computer's complex square 
root subroutine, so that the real part is always posi- 
tive. When the plus sign is used, the wave will be 
called ordinary, and the minus sign will correspond 
to the extraordinary wave. The quantity F, (2.2) 
may be rewritten in the form 



F = Y 2 cos 2 d[(Z 2 -Z 2 ) + (X- 



l) 2 
+ 2iZ{X-\)] 



(2.3) 



where Z c — 



Y sin 2 6 
2 cos 6 ' 



Considering an upgoing wave in the ionosphere, 
initially it will be at a level where X <1, so that the 
quantity F, (2.3), will have a negative imaginary com- 
ponent and will be in the third or fourth quadrant of 
the complex plane. At X=l, F crosses the real axis 
and has the sign of {Z\ — Z 2 ) . Thus if Z c > Z at X = 1 
then VF changes continuously through X=l, and so 
too does the complex refractive index n; but if Z c < Z 
at 1=1 then W jumps discontinuously from the 
negative imaginary to positive imaginary axis and con- 
tinuity of the refractive index n requires reversing 
the upper and lower signs as we cross the X= 1 level. 
Thus an 'initial ordinary' becomes extraordinary, and 
vice versa. 

The complex refractive index n becomes zero when 



X=l-iZ 
for the ordinary wave and 

X=l±Y-iZ 



(2.4) 



(2.5) 



for the extraordinary wave. These reflection points, 
labelled RO, RX+, and RX— respectively occur at com- 
plex heights which are usually well within a wave- 
length of the real axis. 

The ordinary and extraordinary refractive indices 
(and polarizations) become equal at the coupling point 



The four points RO, RX+, RX—, and C are, in general, 
the four branch points of the complex refractive index 
surface. 

In order to locate the branch points in the complex 
plane, it is necessary that X and Z be known analytic 
functions of z = x~\~iy, the complex height. If an 
arbitrary ionospheric model or 'profile' is used, e.g., 
if X and Z are experimentally derived functions of the 
height z, then in order to locate the branch points in 
the complex plane it is necessary in each case to 
approximate the X and Z variations at real heights in 
the neighborhood of these points by analytic functions, 
and then by substitution to solve (2.4) to (2.6) for com- 
plex z. In general the branch point closest to the 
real axis is the one chosen. 

2.2. Location of the Branch Points for Exponential 
Variation of X and Z 

In the present investigation it has been assumed that 
log X and log Z vary linearly in the neighborhood of 
the branch points, so that 



X=X e a * 

Z = Z e- b * 



(2.7) 
(2.8) 



where X and Z are the values at X = 1 or X = 1 ± Y, 
and the constants a and b are determined by measuring 
the slopes of the log X and log Z profiles, or of the 
chords across one or several tabular intervals near 
these points. 

Substituting (2.7) and (2.8) into (2.4) and (2.5) it is 
found that the coordinates (xr, yn) of the reflection 
points, measured from the X=l or X= 1 ± Y levels, 
are given by 



Xr = - log 
a 



cos byR 



cos (a~\-b)yRJ 
1 



log 



[ Zp cos (a 
— Xo si 



i- b) y R 
sin ayR 



(2.9) 



Here yR is necessarily negative and only such values 
of yR are accepted for which the arguments of the sin 
and cos terms are less than 7r/2. Inspection of (2.9) 
shows that if a and b are positive, i.e., if the X profile 
increases and the Z profile decreases with height, then 
xr > 0, i.e., the reflection points will lie at slightly 
greater heights than the X=l or X=lztY levels. 
For a positive and b negative as in the F region where 
electron-ion collisions predominate and consequently 
Z increases with height, x will often be negative. 



2 A second coupling point at 



X = 1 — iZ c — iZ 



which would be located farther away from the real axis, and would hence be less significant 
physically, has been ignored. 
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Substitution of (2.7) and (2.8) into (2.6) yields the 
coordinates (x c , y c ) of the coupling point, measured 
from the X = l level, as given by 



x c - 



log 



cos by c — Z c sin by r 
cos (a + b)y c 



=\\og 



Z cos (a + 6) y c 



Z c cos ay c — sin ay c 



(2.10) 



When Z c > Z at Z= 1, the case of most general interest 
at high frequencies, the positive value of y c is chosen, 
and for Z C <Z at X=l, the negative value of y c . 
In general x c may be positive or negative and x c may be 
greater or less than xro, i.e., the coupling point C 
may be either above or below the ordinary wave re- 
flection point RO. 

Equations (2.9) and (2.10) may be solved numerically 
on a computer. Equating the expressions involving 
y only, allows y to be found by means of a suitable 
zero search procedure, and x is then given explicitly 
as a function of y. 

2.3. Phase Integral Calculation of Absorption 
and Virtual Height 

Having located the main branch points, z , the com- 
plex refractive index must now be integrated from the 
starting level, usually the ground, along suitable con- 
tours up to these points and back again. The complex 
reflection coefficient, measured at the ground, is given 
by 



exp -^ (ct-2 ndz) 



and the total absorption A suffered by the wave is 
(in nepers) 



•■•w 



ndz 



(2.11) 



Representing a pulse or wave packet at the ground, 
by an expression of the type 



JW) exp ff(f -/f^) 



df 



we see that the phase at the ground will be stationary 
with respect to / at a time t = T given by 



2 dr 



f 



ndz 

Jo 



(2.12) 



yielding the virtual height of reflection, h' . 

When X and Z are known analytic functions of z 
it is frequently convenient to choose as contours of 
integration straight lines in the complex plane joining 
the origin to z [Cooper, 1961]. But when X and Z are 
arbitrary tabulated functions of z it is necessary to 
perform the integration mainly along the real axis. 



A convenient path is along the real axis to xo=M(zo) 
and then parallel to the imaginary axis to Zo = Xo + iyoy 
X and Z having the assumed variation in this region 
given by (2.7) and (2.8). This particular choice of 
contours has the advantage of allowing a direct com- 
parison with simple ray theory. 

Broadly speaking, the real axis integration (excluding 
possibly the last wavelength of contour path in the 
reflection region) represents the W.K.B. solution of the 
propagation equations. The 'complex plane integra- 
tion' from the real axis down to the branch point 
Zo, together with perhaps the last wavelength of contour 
path on the real axis, serves as a 'connection formula' 
linking the 'upgoing' to the 'downgoing' W.K.B. 
solutions on the real axis. This connection formula 
represents the asymptotic approximation to the Airy 
integral solution of the propagation equations near the 
branch point where the W.K.B. solutions break down 
[Budden, 1961]. The 'complex plane integration' 
is thus a mathematical procedure for calculating the 
absorption loss and phase change specific to the re- 
flection (or coupling) process where the W.K.B. 
approximations are invalid. 

The four main integration paths are shown sche- 
matically in figure 1. There are two paths correspond- 
ing to the extraordinary waves (X waves) reflected 
at RX + and RX-, respectively, which are important 
above and below the gyrofrequency, respectively. 
The ordinary wave (0 wave) is reflected at RO, and 
the path yielding the Z-trace is that which proceeds 
from ground to the coupling point C in the ordinary 
mode and thence to RX + as extraordinary. We shall 
call the wave traveling along this path the Z-wave. 

Writing the complex refractive index 



f^-iX 



(2.13) 



where /x and \ are tne refractive and absorption indices 
respectively, and substituting Zo = Xo + iyo in (2.11), 
we obtain 



A = 



477/ 



xdx-\ | 

Jo Jxq 



fMdy 



(2.14) 
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Figure 1. 



Key: 



Schematic representation of the four main integration 
paths. 

ordinary waves; extraordinary waves. 
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yo is of course negative and the real axis and complex 
plane absorptions are additive. The Z wave absorp- 
tion may similarly be split into real axis and complex 
plane contributions. 

The virtual height calculation (2.12) is performed by 
calculating 



mf\ ndz =f fidx+ X dy 

L Jo J L Jo Jx 



(2.15) 



for two nearly equal frequencies, thus giving the 
derivative. 



3. Results for Daytime Models 
3.1. Ionospheric Models Used 

In order to demonstrate typical results of the com- 
puter programme two daytime ionospheric models 
have been chosen. In figure 2, curve A shows the 
electron density distribution in the D and E layers 
over Oslo at noon during a sunspot maximum equinox. 
This curve has been proposed by Piggott and Thrane 
[1964]. Curve B shows the variation over Slough at 
noon during a sunspot minimum equinox. The 
Z)-layer curve has been determined by Deeks [1964] 
and the is-layer 'nose' has been tentatively fitted to it 
by Piggott [1964, private communication]. The 
inclination of the earth's field to the vertical at Oslo 
and Slough are = 17° and = 22°, respectively, and 
the gyrofrequencies are/// =1.3 Mc/s and/// =1.28 
Mc/s, respectively. The dashed curve represents the 
mean collision frequency adopted. The two electron 
density distributions contrast cases with relatively 
dense and relatively weak D-region electron densities, 
but also differ significantly in that the slopes, dN/dh 9 
are quite different for most of the reflection levels at 
the frequencies used in our computations. 
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Figure 3. Absorption curves for the Oslo model. 

The X + trace (/<///) is off scale. 
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Figure 4. Virtual height curves for the Oslo model. 









'ogio 17 — * 






3 


4 5 6 7 8 


120 




' 


1 ! 1 1 




- 




J 1 


100 














/ x 


80 


- 




(SLOUGH)/ y^ X 

B*/ yf A N ^ „ 




- 




^ — O (OSLO) N N 


60 




1 


i i i i s 



FIGURE 2. The assumed variation of electron density, N, with height 
at equinox noon. 

Curve A is for Oslo at sunspot maximum and Curve B is for Slough, sunspot minimum. 
Curve C is the assumed height variation of electron collision frequency v. 
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Figure 5. Absorption curves for the Slough model. 
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Figure 6. Virtual height curves for the Slough model. 



3.2. Variation With Frequency of Absorption and 
Virtual Heights 

Figures 3 and 4 show the absorption and virtual 
height variation for the Oslo model, and figures 5 and 
6 repeat these for the Slough model. For all values 
plotted, Z<Z C at X=l. Below the gyrofrequency 
the O- and Z-traces are the least attenuated, and in 
the case of the Oslo model the X + absorption does not 
appear as it is off scale for frequencies above 0.9 Mc/s. 
For both models used, an ionogram might be expected 
to show only an O-trace together with a much fainter 
Z-trace below about 2 Mc/s, while an X-trace would 
only appear near the penetration frequency of the 
O-trace. 

3.3. Complex Plane Contribution to the Total 
Reflection Loss 

The complex plane contribution to the total absorp- 
tion (given by the second term in (2.14)) has been cal- 
culated for all branch points. It never exceeds 
0.05 Np for the X=l + Y reflection in the Slough 
model and 0.015 Np in the Oslo model, and is quite 
negligible compared with the large real height absorp- 
tion in this mode. For Oslo the X= 1 — Y absorption 
in the complex plane (expressed in nepers) is less than 
1 percent of the total for all frequencies above 2 Mc/s, 
and for Slough it constitutes some 5 percent (0.56 Np) 
at 2 Mc/s falling off to about 3 percent (0.18 Np) at 
3 Mc/s. The O-mode complex plane absorption 
nearZ=l is more significant. For Slough it constitutes 
0.48 Np out of a total of 3.37 Np at 0.8 Mc/s, falling 
off to 0.24 Np out of 4.31 at 2.8 Mc/s. The Oslo figures 
are 2.67 out of 14.56 Np at 1 Mc/s falling to 0.26 out of 
3.70 Np at 2.7 Mc/s. The main conclusion is that the 
complex plane absorption is generally significant only 
for the O-wave and especially at the lower frequencies 
which are reflected in regions of higher collision 
frequency. 

3.4. Virtual Height of the Z-Trace 

An interesting feature of the Oslo virtual height 
curves is the appearance of the Z-trace below the 



O-trace at the low-frequency end. This is at first 
sight surprising when analyzed in terms of simple ray 
theory, which considers the Z-mode wave packet to 
travel up to the coupling level near X = 1 in the O-mode 
and then to travel in the A'-mode to the 1 + Y reflection 
level. The first part of the journey might be expected 
to contribute a virtual height equal to that of the 
O-trace, and the second part to add an additional virtual 
height equal to the real height difference between 
the X=l and X = 1 + Y levels plus any retardation in 
this section. The explanation of this anomalous result 
in the light of the phase integral analysis is given in 
section 4. 

3.5. The Z-Trace Intensity as a Function of Magnetic 
Latitude 

The effect of increasing the angle of inclination 
6 between the earth's magnetic field and the vertical 
is to decrease sharply the amplitude of the Z-trace 
by increasing the coupling loss or, in terms of the phase 
integral, by increasing the complex plane contribution 
to the absorption integral, (2.11). We find, for in- 
stance, by using the Slough model distribution for a 
frequency of 1.1 Mc/s, and letting the inclination 6 
take successive values, 11.5°, 17°, 22°, 30°, and 45°, 
that the complex plane contributes an absorption loss 
(for two traversals, upgoing and downgoing) of 0.04, 
1.2, 3.0, 7.1, and 19.2 Np respectively. At 1.5 Mc/s 
the corresponding figures are 0.10, 2.5, 5.5, 12.9, 
and 30.0 Np respectively. 

3.6. Accuracy of the Computations 

The O-mode absorption calculations for the above 
and other ionospheric models have been checked 
against calculations made by a phase integral pro- 
gramme developed by K. G. Budden for the Cambridge 
Edsac II computer. Agreement is good in all cases. 
Slight discrepancies can be traced to the different 
analytic approximation made to the electron distri- 
bution near branch points. In our programme the 
curves are assumed to be exponential over short 
distances; in the Edsac II programme various other 
methods of curve fitting were used. This has an effect 
both on the computed location of branch points, and 
also on the interpolated values of X and Z between 
relatively widely spaced tabular intervals. The 
Edsac II programme checks against Budden's and 
against Pitteway's full wave computations. 

Virtual heights for the above and other models for 
the ordinary and extraordinary modes have been com- 
pared with results given by a Mercury computer 
programme developed by E. B. Thrane in Oslo. 
Agreement is again good in all cases. 

4. Calculation of Virtual Heights by Simple 
Ray Theory and by the Phase-Integral 
Method 

4.1. Phase-Integral Correction to Ray Theory 
In section 2 it was shown (2.12) that the virtual height 
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of reflection is given by the expression 



h' =^M 
df 



f\ n (z) dz\ 



(4.1) 



Since the upper limit of integration z is a function of 
/this may be written 



h'=% 



where 



n'dz+fn{zo) —£ 
Jo d f. 



n'-^W 



(4.2) 



(4.3) 



defining the complex group refractive index. 

Since z represents a reflection branch point, n(zo) 
= 0, and hence we obtain 



h' = % n'dz. 

Jo 



(4.4) 



Equations (4.1) and (4.4) are two equivalent phase 
integral expressions for the virtual height, but the 
method of evaluation will be different for each. 
Consider (4.4) first. Since n and hence n' are known 
analytic functions of frequency, n' may be calculated 
for all points on the integration path and then inte- 
grated to zq. This method is adopted by Cooper 
[1961]. If we compare this method with the ray theory 
formulation 



h'=\ fi dx 



(4.5) 



(where /jl 1 ' =-&(n') is the group refractive index) 
and ignore for a moment the fact that Xo in (4.5) is 
not quite equal to ffl (z ) of (4.4), we find one form of the 
phase integral correction A/i', to the ray theory esti- 
mate. If OAiRi (fig. 7) is the integration path chosen, 
then 



M'=^f n'dz=\ y x'dy (4.6) 

JA l R 1 AxRiJo 

where x' is the group absorption index. The evalua- 
tion of this correction term necessitates doing a com- 
plex plane integration and is not readily evaluated 
if a phase integral computer programme is not avail- 
able. 

The alternative method of evaluating h' as adopted 
in the present investigation is based on (4.1) and gives 
a clearer insight into the cause of retardation, yielding 
also a more useful expression for the correction to 
ray theory. The quantity 



n 



. Jo 



n{z)dz 
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X (fi) X (f 2 ) 
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Figure 7. Contours in virtual height calculation. 

OAiRi represents an integration path to reflection point Ri for a frequency /i and OA2R2 
is the same for a slightly higher frequency fo. 



is evaluated for two neighboring frequencies, f\ and 
/2, the integration paths being represented schemat- 
ically by OA\R\ and OA 2 Rz, respectively (fig. 7). As 
fi—*fi the virtual height is given by 



h'=m 



f 2 I n dz—fi n dz\ l(f*—fi 



{fo(f*)-fili(fi)}dx+fi\ fx(f 2 )dx 

+ V 2 [ X (f*)4r 

J A 2 R 2 

-/. f x(fi)dy) 

J A X R X 



Hfz-fi) .... (4.7) 



Now the last bracketed term in the numerator repre- 
senting the difference between the complex plane 
integrations at neighboring frequencies is usually 
negligible compared with the first two terms. For 
instance for the Slough model in the entire frequency 
range considered, this contribution never exceeds 
0.2 km for the 0, X, or Z waves (unless the absorption 
is in excess of 100 Np). Hence (4.7) in the limit re- 
duces to 

h'=J °n'dx+f(JL-jjr- (4.8) 

This equation may have been derived directly from 
(4.2) if the complex plane contribution were ignored. 
The second term on the right-hand side of (4.8) 



Ah' = fr^ 



(4.9) 



which we shall call the 'real axis correction term' 
represents quite closely the difference between the 
phase integral and the ray method estimates of vir- 
tual height. We have thus reduced the complex 
plane correction to ray theory to a term involving 
quantities on the real axis only, the real axis coordi- 
nate Xo of the branch point and its frequency variation 
being all that is required. 
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4.2. Magnitude of the 'Real Axis Correction' to Ray 
Theory 

Referring to (4.9) it is seen that Ah' is zero only 
when there are no collisions and /jl(xo) = 0. Its mag- 
nitude will be enhanced when dN/dh is small, i.e., 
near a peak or in a 'ledge' of the electron density 
distribution. For the E region ix(x ) is generally 
appreciably larger for the O-wave (i.e., near X=l) 
than for the X waves (i.e., near X= l±Y) except when 
the latter are heavily attenuated, so that the correction 
term is usually more important for the O-mode reflec- 
tion. Neglect of this term in a true height analysis 
will cause overestimation of layer heights. 

The term dx ldf and Ah' have been calculated for 
a number of frequencies for the Oslo model, and values 
of Ah' for the O-mode are shown in table 1. Values 
of /n(xo) are also given for the X=l and X— \±Y 
reflection points, labeled RO, RX+, and RX-, re- 
spectively. The quantity Ax R is the amount that x 
lies beyond the X=l level. 



Table 1. Variation of the 'real axis correction term' Ah' with 
frequency for the Oslo equinox noon value 

(Symbols are defined in the text) 



/ 


fi(x ) 


Ax« 


AA' 


Mcls 


RX- 


RO 


KX + 


(O-mode) 


1.0 
1.2 
2.0 
2.5 
3.0 




0.41 
.30 
.20 
.16 
.11 


0.056 
.053 
.041 
.030 
.020 


km 

0.0066 
.029 
.022 
.001 
.007 


km 

7.75 




2.06 


0.11 
.076 
.060 


0.48 

.92 

2.65 



We observe that Ah' may be quite appreciable in 
the E layer. It is seen that x is usually very close to 
the ray theory reflection level X=l (the same is true 
for the X= 1 ±Y levels). This means that little error 
is introduced by performing the real axis integration 
to the X=l and X=l±Y levels. It also means that 

the quantity —^ appearing in the formula for Ah' can 

be approximated by calculating the frequency varia- 
tion of the X= 1 (or X= 1± Y) levels, which is a rela- 
tively simple procedure. 

We note finally that near the peak of a noon F 2 layer, 
Ah' is of the order of 5 km, but lower down it is very 
much smaller. 



4.3. Virtual Height of the Z-Trace 

The virtual height of the Z-trace may be calculated 
by similar considerations. We calculate the expression 
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FIGURE 8. Integration contours for determining the virtual height 
of the Z-trace. 
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for two neighboring frequencies and hence for two 
slightly different integration paths (fig. 8), and then 
calculate h' as in (4.1). This yields three main terms 
as in (4.7), the first corresponding to a simple ray 
theory calculation, the second representing a 'real 
axis correction term', and the third representing the 
difference between two pairs of complex plane inte- 
grations, which again turns out to be negligible com- 
pared with the other two terms. 

The 'real axis correction term' turns out to be 



Ah'- 



/JL (C) — fJL x C) 



— +flx (R)f w (4 . 10 ) 



where the subscripts o and x refer to the ordinary and 
extraordinary polarizations, and C and R refer to the 
coupling and RX-\- reflection points respectively. 
The first term 



Mo(C) 



»AC)]ff 



is the real axis correction term corresponding to the 
coupling point integration. It is made up of two parts, 
the first 



Ho(C)f 



dx c 

w 



being roughly equal to the usual O-wave real axis 
correction, and the second 

-MC)/f 

whose effect is to decrease the Z-trace virtual height by 
this amount with respect to the O-trace. In the Oslo 
model at 1000 kc/s where the O-trace is 11 km above 
the Z-trace, this term is equal to —20.5 km and is 
only partially compensated by the additional integra- 
tion to the RX-\- reflection point. 

In physical terms one might say that the fact that a 
Z-trace may appear below an O-trace expresses the 
fact that the concept of a wave packet traveling to a 
fixed point near the X= 1 level and there being either 
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reflected, or by coupling continuing to the RX+ level, 
has only limited validity. For, in fact, each frequency 
component of the wave packet travels along a slightly 
different path and undergoes coupling or reflection at 
slightly different levels. This gives rise to a correction 
term in the virtual height determination which can 
produce, for instance, the appearance of a Z-trace 
below an (9-trace. 



5. Z-Trace Both Above and Below the 
Gyrofrequency 

In the F-region Z c > Z, unless 6 ~ 0, and the com- 
plex plane absorption loss for the Z-trace is usually 
very large. Using typical F-region distributions of 
electron density and collision frequency [Johnson, 
1961] values of complex plane absorption above 100 
Np are common, even for 8= 17°. Hence the appear- 
ance of F-layer Z-traces in high latitude ionograms 
cannot be interpreted by means of simple magneto- 
ionic coupling in a horizontally stratified ionosphere 
as described here. Budden [1961] comes to the same 
conclusion. 

We have seen that the F-region Z-trace, which is 
normally ignored in theoretical analyses for frequencies 
below the gyrofrequency, in fact continues smoothly 
downwards to frequencies well below the gyrofre- 
quency where it becomes one of the least attenuated 
traces. For the daytime E layer when /</// the 
Z-trace coupling region lies well below the F-layer 
peak and Z is comparable with Z c . As we proceed to 
lower frequencies the RO and C branch points appear 
at lower heights, where Z is larger, and the coupling 
point C (fig. 1) moves down until it reaches the real 
axis when Z = Z C at X=l. For the Slough model this 
would occur at 0.25 Mc/s, and for Oslo at 0.8 Mc/s. 
At this stage the upgoing ordinary becomes an upgoing 
extraordinary, with the coupling loss having just be- 
come zero. As /is decreased further, Z c < Z atX= 1, 
and the upgoing ordinary becomes an upgoing extraor- 
dinary at X = 1 without coupling, and after reflection 
eventually reaches ground with the ordinary wave 
polarization. The 'initial Ordinary' trace thus pro- 
duced would hence be a continuation of the Z-trace 
and would usually be the strongest trace in this fre- 
quency region. Piggott, Pitteway, and Thrane [1964] 
have indeed found in their full-wave analysis at low 
frequencies that very little reflection takes place at 
any height besides the X=l + Y level. A second, 
weaker trace produced by reflection at the X=\+Y level 
would be due to the (initial) extraordinary wave, which 
by coupling would continue past the X=l level as 
extraordinary, involving a discontinuity, i.e., the phase 
intregal contour path would have to go below the 
coupling point in the negative imaginary plane [Cooper, 
1961]. 

Figures 9 and 10 are daytime ionograms taken at 
Fort Churchill (0 = 6°) showing O- and Z-traces for 
the E and F regions. At the low frequency end of 
figure 9, we see that the Z-trace is the only remaining 
trace. In figure 10 we see that for a small frequency 
range the Z-trace is unmistakably below the O-trace. 
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Figure 9. Ionogram taken at Fort Churchill at 5:30 p.m. on 
Sept. 13, 1958. 

At the extreme low frequency end only the Z-trace is present. 
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FIGURE 10. Ionogram taken at Fort Churchill at 10:01 a.m. on 
Sept. 6, 1958. 

The Z-trace is below the O-trace for a small frequency range. 



6. Conclusions 

The use of a phase integral computer programme to 
produce absorption and virtual height versus frequency 
curves allows direct comparison with ionograms and 
pulse absorption measurements. We have shown 
that the Z-trace, which is normally ignored in the the- 
oretical analysis for frequencies below the gyrofre- 
quency, continues smoothly, in fact, to frequencies 
well below the gyrofrequency, where it becomes one 
of the two least attenuated traces. 

It has been shown that the phase integral calcu- 
lation of virtual heights can be approximated quite 
accurately by using the simple ray theory calculation 
together with a 'real axis correction term' only (4.9), 
requiring only the location of the branch point in the 
complex plane at two neighboring frequencies. In 
many cases this correction term can be approximated 
by calculating the simple ray theory reflection heights 
where X=l and X=1±Y, for two neighboring fre- 
quencies, together with the value of \x at these heights. 

Finally it has been shown that the Z-trace virtual 
height calculation will include a term 



518 



which in some cases may be sufficient to produce a 
Z-trace at a lower virtual height than that of the 
O-trace. 

The 'real axis correction term' is usually more sig- 
nificant in the E region than in the F region, where it 
is doubtful whether it would be an important factor 
in explaining the overestimation of heights in electron 
density profiles which has sometimes led to a crossing 
over of 'topside' and 'bottomside' F-region profiles. 



This communication is published by kind permis- 
sion of the Director of the D.S.I.R. Radio Research 
Station. The ionograms for Fort Churchill are pub- 
lished by permission of the World Data Centre for the 
IGY in Slough. The author is, indebted to Mr. W. R. 
Piggott for valuable discussions in the later stages 
of this work, and in particular for pointing out the 
existence of the £"-region Z-traces for high latitude 
stations. 
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